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I.  INTRODUCTION 

In  recent  years,  a  number  of  authors  have  considered  the  problem 
of  wake  collapse  in  an  incompressible,  stratified  fluid.  These  treat¬ 
ments  nave  usually  involved  a  version  of  the  linear  Boussinesq  approx¬ 
imation,  though  some  numerical  work  has  appeared.  While  there  seems 
to  be  a  consensus  about  the  phenomenology  in  the  case  of  fully  mixed 
wake  (for  which  the  linear  approximation  is,  of  course,  invalid),  it 
has  seemed  to  us  useful  to  treat  a  particularly  simple  version  of  the 
linear  problem,  for  which  an  exact  solution  is  obtainable.  The  re¬ 
sults  contain  some  surprises,  and  provide  some  insight  into  the  limi¬ 
tations  of  the  linear  treatment. 
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II.  BASIC  EQUATIONS 


We  study  a  linearly  stratified  fluid,  with  no  boundaries,  whose 
unperturbed  density  PQ(z)  in  the  vicinity  of  the  region  of  interest 
varies  only  in  the  vertical  z -direction,  and  gradually  enough  to 
justify  a  linear  approximation.  Thus 

°o<z)  ~  °oo  ■  BZ  (1) 


where  poo  is  a  "mean”  density  for  the  problem.  The  initial  value 
problem  is  obtained  by  perturbing  the  density  slope  inside  a  cylinder 
of  radius  a 


6p(t=0)  =  ez  r  <  a 

=  0  r  >  a 


(2) 


where  the  fully  mixed  case  can  be  obtained  formally  by  setting  e  =  8. 
This  state  of  affairs  is  illustrated  in  Fig.  1.  We  are  working  in 
a  two-dimensional  coordinate  system  with  axes  z  (vertical)  and  x 
(horizontal),  with  everything  assumed  independent  of  y.  We  will 
also  use  polar  coordinates  (r,a)  in  the  plane  of  the  problem,  with 
the  polar  axis  vertical. 

Newton's  second  law,  linearized,  is 

Sv 

P  TT  =  -  7P  -  pg  S  (3) 

OL  **  r* 


where  v  is  the  fluid  velocity  (in  the  x-z  plane),  p  is  the  pressure, 
g  the  acceleration  of  gravity,  and  z  a  vertical  unit  vector.  Since 
the  fluid  is  assumed  incompressible,  we  have 
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divv-^r+v‘7p-0  (4) 

which  allows  us  to  introduce  the  stream  function  ^ ,  a  vector  in  the 
y-direction,  defined  by 


v  -  curl  ^ 


(5) 


The  curl  of  Eq.  3  yields 


poo  kt (curl  curl  V  = 


g  7p  x 


(6) 


where  the  insertion  of  pQo  is  usually  called  the  Boussinesq  approxi¬ 
mation.  It  amounts  to  taking  into  account  any  density  gradient  that 
leads  to  a  force  (is  multiplied  by  g),  but  no  other.  Differentiating 
Eq.  6  with  respect  to  time  leads,  in  view  of  Eq.  4,  to 


-  N 


2 


1? 


(7) 


where  we  have  used  the  linear  approximation,  and  have  defined  the 
Vais ala  frequency 


o  dp 

N2  ,  -  —2  „  gB/p 

pQ  dz  w  yK/Hoo 

which  we  will  treat  as  a  constant  hereafter.  (Obviously,  we  could 
have  chosen  pQ  exponential,  instead  of  as  given  in  Eq.  1.) 

Clearly,  in  view  of  Eq.  4,  ftp  can  be  obtained  from  *  ,  so  that 
Eq.  7,  with  the  initial  condition  Eq.  2,  is  our  problem. 

It  begs  for  Fourier  analysis,  and  the  (exact)  answer  is 
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where  the  ±  means  that  the  integral  should  be  evaluated  with  each  of 
the  two  signs,  and  *-he  results  averaged.  The  corresponding  expression 
for  jk  is  obtained  from  Eq.  9  by  multiplying  the  integrand  by  2N/k|S, 
and  taking  the  difference  of  integrals  for  the  two  signs.  (There  are 
other  simpler  ways  to  get  the  velocity  distribution.) 

It  is  worth  reemphasizing  that  Eq.  9  is  the  exact  solution  to  a 
linear  problem. 


III.  QUADRATURES 


The  evaluation  of  Eq.  9  is  greatly  aided  by  the  observation  that 

the  expression  in  the  integrand  involving  the  square  roots  admits  of 

an  expansion  in  terms  or  Tschebysheff  polynomials  of  the  second  kind. 
2  2  2 

Explicitly,  if  w  =  u  +  v  -  2uv  cos  y. 
then 


w 


os  v) 

m  =  0 


v  m 


(10) 


where 


A  (cos  Y>  =  *ln  l^hl 
m  T  sm  Y 


(11) 


Thus,  after  a  bit  of  arithmetic,  we  find 


6o  = 


2za‘ 


^.,i(Nt) 


^^(-1^  (24+1)  C°S’C(P^1)°1 


i=0 


/' 


J2(ka)  J24+l(kr)  dk  (12) 


Before  proceeding  to  the  final  form,  some  observations  about  Eq.  12 
are  in  order. 
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First  Aid  foremost,  inside  the  original  cylinder,  where  r  <  a, 
only  the  fir.it  term  in  Eq.  12  is  non-vanishing.  Thus,  exactly. 

vMNt) 

6p  -  2ze  — — ;  r  <  a  (13) 

The  original  linear  (proportional  to  z)  perturbation  inside  the 
cylinder  remains  linear,  overshoots  its  correction,  anJ  finally  damps 
out. 

The  fluid  velocities  for  r  <  a  are  most  easily  obtained  from 

Eq.  13,  by  noting  that  v  =  —  ^*(6p),  so  that 

z  p  ot 


vz  =  -  J2(Nt)  (14a) 

In  view  of  the  incompressibility  condition  (Eq.  4),  we  also  have 

\  =  221  Js(KC)  (14b) 

so  that  the  fluid  particles  inside  the  original  cylinder  move  on  right 
hyperbolae,  overshooting  their  ultimate  positions  on  the  first  pass. 
The  ultimate  displacement  of  a  fluid  particle  that  starts  inside  the 
original  cylinder  is  given  by 


A*  = 


£21 

B 


A  z  =  - 


*Z 

8 


os) 


so  that  the  fluid  particle  moves  just  far  enough  down  its  own  hyper¬ 
bola  to  reach  the  level  which  is  appropriate  for  it.  In  particular, 
the  original  circular  cross  section  of  the  cylinder  deforms  into  an 

ellipse  of  semi-major  axis  a(l  +  ■*■)  and  semi-minor  axis  a(l  -  -f-), 

P  p 

after  oscillating  around  this  shape. 
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We  now  turn  to  the  behavior  of  the  fluid  for  r  >  a,  for  which 
all  terms  except  the  first  in  Eq.  12  contribute.  We  need,  for  r  >  a, 
and  for  t  2  1, 


/ 


J2(ka)JH+i(kr)  dk  =  — j 


4  a 


[p<..i'2pt + + \  -  H  <io> 


where  the  are  conventional  Legendre  polynomials,  whose  argument  is 

1  -  The  expression  in  square  brackets  can  be  written  in  other 

r^ 

forms,  using  the  many  identities  involving  the  Legendre  polynomials. 
If  we  define 


*,_!<?>  -  P.+1(S> 

G*<5>  H  Vl(^  -  2Pt<5>  +  +  ^  h h  T A  . 


2  (i  -  ?)  P.(s) 


then  the  first  two  are  given  by 


pt-i<s>  -  'W?>‘ 
1 


21  + 


4- 7) 


2a2  \  _  4a4 


2  /  4 

r  /  r 


12a4  16a6 


,  etc. 


We  have,  then, 

6p 


(17) 


(18a) 


(18b) 


= !*E  J2tNt(Nt>  Gt(i  -  *4) 

1=1  \  r  / 

! -  E  c05i°s^1?aJ  ^f1  •  *4) 

1=1 


2z 

r 


(19) 

(20) 
r  >  a 
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where  (again,  with  5  1 - ^  ' 

r 

H.(?)  =  -  G  (?) 

c  4a'  1 

-  (24*i)Pt(?>  ♦£r(*W?)  •  Vi(o)  (21) 

It  is  now  interesting  to  observe  that,  as  r  -  a,  H^(?)  -  (-])^(24 
so  that 


6p(r-«a  )  =  -2zc 


cos[(24+l)o]  J24+l(Nt) 
cos  a  Nt 


(24+1) 


(22) 


4=1 


Comparing  this  with  Eq.  13,  we  find  an  exact  expression  for  the 
density  discontinuity  across  the  surface  of  the  cylinder 


A 0  =  6p(  r-a")  -  6p(r-a  +  ) 


(Nt) 


(23) 


4=0 


This  sum  can  be  extended  to  -®,  arid  carried  out  exactly. 

We  find 

AP  =  ez  cos(Nt  sin  or)  (24) 

an  unexpected  result.  From  Eqs.  24  and  13,  we  find 

+  r  jx(Nt)  i 

6p(r-a  )  =  cz  2  -=rjr - cos(Nt  sin  <*) 


(25) 


which  does  not  go  to  zero  as  t  - ■*>  ,  but  which  oscillates  more  and  more 
rapidly  as  a  function  of  a*  Such  behavior  is  indicative  of  an  insta¬ 
bility  in  a  more  realistic  calculation,  thouah  the  relation  to  the 
expected  Rayleigh  instability  on  the  surface  is  unclear. 

Wo  turn  finally  to  the  behavior  at  large  distances,  for  which 
we  need  the  behavior  of  the  H.  for  5  near  unity 


Ht<!> 


-  ( 24+1)  J, 


■  “ 

(*W)  fj  i 


r  »  a 


(26) 


in  which  we  have  also  assumed  1  «  i  «  —r  for  convenience.  It  will 

a* 

be  seen  below  that  this  is  the  range  of  interest  for  4.  This  leads 
to 


6p 


<-ir(2t-rl) 


cos  Itailiial  'hiik 
cos  a  Nt 


(Nt) 


jwi)  f] 


t-l 


(27) 


(Recall,  as  always,  that  v  can  be  obtained  from  6P  through  v  =  —  ( 6 P ) , 

Z  Z  p  QL 

and  v  from  v  through  the  incompressibility  condition  Eq.  4.) 

X  z 

For  any  given  point  (r,a)  and  time,  the  terms  in  the  sum  (Eq.  27) 
increase  in  magnitude  with  l,  and  oscillate  rapidly  in  sign,  so  it  is 
appropriate  to  look  for  a  "constant  phase"  value  of  t  as  a  means  of 
estimating  the  sum.  By  either  this  method,  or  by  a  direct  saddle- 
point  integration  in  Eq.  9,  we  find  that 

6P  -  ---J-  tan  a  sin(Nt  cos  a)  J?  sin  (28) 


for  r  »  a.  This  represents  a  pulse  of  frequency  N  cos  a,  which 


passes  a  point  (r,  a)  at  a  time  t  ~ 


The  denominator  is 


Na  sin  a 

the  group  velocity  of  a  wave  of  wave-number  ~  1/a,  in  the  correct 
direction,  where  N  cos  0  is  the  frequency  of  such  a  wave.  (Recall 


that  the  group  and  phase  velocities  of  Vais ala  waves  are  mutually 
orthogonal. ) 
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IV.  CONCLUSIONS  AND  CAVEATS 


The  behavior  of  the  solution  (Eq.  28)  for  large  distances  from 
the  source  contains  no  surprises,  and  represents  the  radiation  from 
the  source  of  the  expected  pulse  of  Vais ala  waves,  necessary  to  get 
rid  of  the  energy  stored  in  the  initial  perturbation. 

On  the  other  hand,  the  behavior  near  the  perturbation  has  a 
number  of  unrealistic  features.  The  rapidly  oscillating  (in  space) 
behavior  just  outside  the  original  cylinder  at  long  times  cannot 
appear  in  the  solution  of  a  realistic  hydrodynamic  problem.  Since 
the  solution  is  exact,  this  difficulty  must  be  ascribed  to  the  model. 
In  addition,  the  velocity  of  the  fluid  displays  an  infinite  shear  on 
the  surface  of  the  cylinder,  and  this  too,  in  a  realistic  problem, 
would  lead  to  Rayleigh  instability.  It  would  seem  foolhardy,  there¬ 
fore,  to  set  c  =  g,  and  use  a  linear  treatment  to  study  the  fully 
mixed  wake  collapse  problem. 
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